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Abstract 

We discuss a matrix public key cryptosystem and its numerical implementation. 


1 Introduction 

In a symmetrical key cryptosystem, such as AES (Advanced Encryption Standard), two users Alice 
and Bob must first agree on a common secret key [1]. If Alice communicates the secret key to 
Bob, a third party. Eve, might intercept the key, and decrypt the messages. In order to avoid such 
a situation one can use an asymmetric public key cryptosystem, which provides a mechanism to 
securely exchange information via open networks [1]. 

In public key cryptography a user has a pair of cryptographic keys, consisting on a widely 
distributed public key and a secret private key. These keys are related through a hard mathematical 
inversion problem, such that the private key cannot be practically derived from the public key. 
The two main directions of public key cryptography are the public key encryption and the digital 
signatures. Public key encryption is used to ensure confidentiality. In this case, the data encrypted 
with the public key can only be decrypted with the corresponding private key. Digital signatures 
are used to ensure authenticity or prevent repudiation. In this case, any message signed with a 
user’s private key can be verified by anyone who has access to the user’s public key, proving the 
authenticity of the message. 

A standard implementation of public key cryptography is based on the Diffie-Hellman (DH) key 
agreement protocol [2]. The protocol allows two users to exchange a secret key over an insecure 
communication channel. The platform of the DH protocol is the multiplicative group Zp of integers 
modulo a prime p. The DH protocol can be described as following: 

1. Alice and Bob agree upon the public integer g G Zp. 

2. Alice chooses the secret integer a. 

3. Alice computes A = g°' modp, and publishes A. 

4. Bob chooses the secret integer h. 

5. Bob computes B = g^modp, and publishes B. 

6. Alice computes the secret integer Ka = B°^ modp = g^°' modp. 

7. Bob computes the secret integer Kb = A^modp = g°'^m.odp. 
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It is obvious that both Alice and Bob calculate the same integer K = Ka = Kb, which then can 
be used as a secret shared key for symmetric encryption. 

Assuming that the eavesdropper Eve knows p,g,A and B, she needs to compute the secret key 
K, that is to solve the discrete logarithm problem: 

A = g'^m.odp, (1) 

for the unknown a. If p is a very large prime of at least 300 digits, and a and b are at least 100 
digits long, then the problem becomes computationally hard (exponential time in log p), and it 
is considered infeasible. For maximum security it is also recommended that p is a safeprime, i.e. 
{p — l)/2 is also a prime, and g is a primitive root of p. This protocol is secure, assuming that it is 
infeasible to factor a large integer composed of two or more large prime factors. 

Here, we discuss a public key cryptosystem, which extends the key exchange to a matrix proto¬ 
col, and avoids the traditional number theory approach based on the discrete logarithm problem. 
The key exchange mechanism is based on a commutative function defined as the product of two 
matrix polynomials. We provide a detailed description of this encoding mechanism, and a practical 
numerical implementation. 


2 Commutativity of Matrix Polynomials 

Let us consider the complex matrix X G , and two complex matrix polynomials of the form: 


and 


P{X, a) : X ^ 

(2) 

M 


PiX,a) = 

(3) 

m=l 


P(X, b) : xC^ ^ 

(4) 

K 



P{X,b) = '£bkX^, 


k=l 


where a G and b G C^. 

One can easily show that the product of these polynomials is commutative: 


M 


K 


\k=l 


P{X,a)P{X,b) = Y. ^raX^ 1 ( Y.^kX'^ 

\m=l 
M K 


K 


m=l k=l 
, / M 


= Y ^kX^ E 

Vfc=l / Vm=l / 

However, because the matrix product is non-commutative, we also have: 

P{X,a)P{Y,b)^P{Y,b)P{X,a), 
if X ^ y G and a £ , b £ . 


(5) 


( 6 ) 


(7) 
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3 Key Exchange Protocol 

Using the above property, the proposed key exchange protocol can be formulated as following: 

1. Alice chooses the secret vectors a € and d G (Alice’s private key). 

2. Alice randomly generates and publishes the matrix U € (Alice’s matrix public key). 

3. Bob chooses the secret vectors b G and b G (Bob’s private key). 

4. Bob randomly generates and publishes the matrix V G (Bob’s matrix public key). 

5. Alice computes and publishes the matrix A = P{U,a)P{V,d) (Alice’s public key). 

6. Bob computes and publishes the matrix B = P{U,b)P{V,b) (Bob’s public key). 

7. Alice calculates the secret matrix Sa = P(U,a)BP{V,d) (Alice’s secret key). 

8. Bob calculates the secret matrix Sh = P(U,b)AP(y,b) (Bob’s secret key). 

One can see that both Alice and Bob obtain the same secret key S = Sa = Sb, since the matrix 
polynomials satisfy the commutativity property. Also we assume that we always have U ^ V. 

Assuming that the eavesdropper Eve knows U, V and A, B, the hard problem is to compute S, 
which requires the unknown private key vectors a, d or b, b. We should note that the number of 
elements (Mi, M 2 , Ji, J 2 ) in these vectors is also unknown. 

One may attempt to find these quantities directly from A and B. For example we have: 

Ml M2 K 

^rn~anU^V^ = afc'I/W. (8) 

m=ln=l k=l 

Thus, we obtain a system of equations,T q; = A, with the unknown ak = am^n, k = {m,n). The 
columns of the matrix T, and the matrix A are written as vectors by stacking 

their constituent columns. In total there are K = Mi M 2 columns in T, and each column has N'^ 
elements. Therefore, the resulted system is badly underdermined if Mi, M 2 ^ N. Thus, we should 
consider that Mi, M 2 , Ji, J 2 ^ in order to increase the security. 

4 Practical Implementation Aspects 

The goal of this section is to discuss some practical aspects regarding the numerical implementation 
of the proposed matrix public key cryptosystem. We choose the C language for implementation, and 
the standard GCC compiler on a desktop running Ubuntu 14.04 LTS. We focus only on building a 
simulation environment, where we can test the functionality of the algorithms. The code is listed 
in the Appendix 1, and here we provide a detailed discussion of its main components. The required 
compilation and run steps are: 

gcc -03 -Im keys.c -o keys 
./keys > results.txt 

The program starts by setting the size of the required square matrices. We choose the size value 
N = 4, which means that all the matrices will contain NN = N*N = 16 double precision complex 
numbers. This value of N is enough to obtain a secret key with the length K = NN*8 = 128 bytes. 
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Since this is a simulation only, we choose also to use the standard random number generator from 
C, randO, in order to initialize the variables in the problem. Of course in a real application rand() 
is a poor choice, and should be avoided, however our goal here is to simply show that the proposed 
algorithm works. For a secure implementation one can use successive applications of secure hash 
functions (sha256), or secure ciphers (AES), to generate the required pseudo-random numbers. 

The private keys, a, d, b, b, are stored in the byte arrays a, aa, b, bb, and they are generated 
with the private_key(. . .) function. The private keys are stored as byte arrays of length Ml*16, 
M2*16, Jl*16, J2*16. Longer private keys will provide a better security. For simulation purposes, 
the Ml, M2, Jl, J2 values are randomly generated in the range [NN, 2*NN-1], such that the condition 
Ml , M2 , Jl, J 2 ^ N is marginally satisfied (one can make them much longer if necessary). 

Each double precision complex number needs 16 bytes to store it. The real and imaginary 
part of the complex numbers are generated with the drandO function. This function is designed 
to generate “dense” double precision random numbers in the ( 0 , 1 ) interval. For each decimal in 
the mantissa, we use a separate rand() call, until the precision gets below 2 • 10 “^®, which is the 
machine epsilon for double precision numbers. 

The matrix public keys, U and V, are stored in the byte arrays U and V, and are generated with 
the matrix_public_key (. . .) function. Each array holds NN double precision complex numbers in 
NN*16 bytes, using the same approach described above for the private_key (. . .) function. 

The public keys, A and B, are stored in the byte arrays A and B, and are generated with the 
public_key(. . .) function. Since these are also double precision complex matrices with the size 
NN, they will need NN*16 bytes of storage. 

The power matrices are also normalized as following: 

(9) 

where ||.||^ is the Frobenius norm. Thus, all the polynomials are rewritten as: 

M 

P{U,a) = Y,am\\U^\\-p^U^, ( 10 ) 

m=l 

This way, the elements of the power matrices will be in the same “range”, while the polynomial 
commuting properties described before are preserved. Also, since the order U, V or V, U cannot 
be specified exactly, we made the convention that if the first element of U is smaller than the first 
element of V, then the order U, V is switched to V, U. This makes irrelevant the order in which 
the arguments, U and V, are passed to the functions. 

The secret keys, Sa and Sh, are stored in the byte arrays Sa and Sb with a length K = NN*4, 
and they are computed with the secret_key(. . .) function. The secret keys are extracted from 
the matrices Sa and Sh, each containing NN elements. The computation of these matrices is affected 
by rounding errors, due to the finite floating point number representation, and inherently will lead 
to Sa 7 ^ Sh, which of course is undesirable. In order to counter the accumulation of the floating 
point rounding errors, we extract the significand of each double precision number, and from the 
significand we extract a 4 byte unsigned integer (see the code for more details). Thus, from each 
double precision complex number we extract two unsigned integers in the range [0,2^^ — 1]. The 
128 bytes corresponding to these 32 integers form the secret keys, and are stored into the byte 
arrays Sa and Sb. 

For illustration purposes, the results obtained for one instance run are given in Appendix 2. One 
can see that both Alice and Bob compute the same secret key (Sa = Sb), using completely different 
secret and public keys. 
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5 Conclusion 


In conclusion, we have presented a matrix public key cryptosystem, which avoids the cumbersome 
key generation of the traditional number theory approach. The key exchange mechanism is ensured 
by the commutative property of the product of two matrix polynomials. Also, we have provided a 
detailed description of this encoding mechanism, and a practical numerical implementation. 
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Appendix 1 

Below we give the C code of the proposed key exchange protocol. 

// 

// keys.c, C implementation of the key exchange protocol 
// (c) 2015, M. Andrecut 
// 

#include <stdlib.h> 

#include <stdio.h> 

#include <math.h> 

#include <complex.h> 

#include <time.h> 

double drandOf 
double d = 0.0; 
double s = 1.0; 
do { 

s /= RAND_MAX + 1.0; 
d += randO * s; 

}while(s > 2E-16); 
return d; 

} 

void cprod(unsigned N, double complex *X, 
double complex *Y, double complex *R){ 
unsigned i, j, n, NN = N*N; 
for(j=0; j<N; j++){ 
for(i=0; i<N; i++){ 

R[j*N + i] = 0; 
for(n=0; n<N; n++){ 

R[j*N + i] += X [n*N + i] *Y [j *N + n] ; 

} 

} 

} 

} 
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double fnorm(unsigned N, double complex *x){ 
double norm = 0, d; 
unsigned n; 
for(n=0; n<N; n++){ 
d = cabs(x[n]) ; 
norm += d*d; 

} 

return sqrt(norm); 

} 

void secret_key(unsigned Ml, unsigned M2, 
unsigned N, unsigned K, 
unsigned char *a, unsigned char *aa, 
unsigned char *U, unsigned char *V, 
unsigned char *B, unsigned char *S){ 
unsigned m, n, i, j, NN = N*N; 

double complex *u = calloc(NN, sizeof(double complex)); 
double complex *v = calloc(NN, sizeof(double complex)); 
double complex *pu = malloc(NN*sizeof(double complex)); 
double complex *pv = malloc(NN*sizeof(double complex)); 
if (U[0] >= V [0] ){ 

for(n=0; n<NN; n++){ 

unsigned char* x = malloc(16); 
for(m=0; m<16; m++) x[m] = U[n*16+m]; 
pu[n] = u[n] = *(double complex*)x; 

} 

for(n=0; n<NN; n++){ 

unsigned char* x = malloc(16); 
for(m=0; m<16; m++) x[m] = V[n*16+m]; 
pv[n] = V[n] = *(double complex*)x; 

} 

} 

else{ 

for(n=0; n<NN; n++){ 

unsigned char* x = malloc(16); 
for(m=0; m<16; m++) x[m] = V[n*16+m]; 
pu[n] = u[n] = *(double complex*)x; 

} 

for(n=0; n<NN; n++){ 

unsigned char* x = malloc(16); 
for(m=0; m<16; m++) x[m] = U[n*16+m]; 
pv[n] = V[n] = *(double complex*)x; 

} 

} 

double complex *c = calloc(Ml, sizeof(double complex)); 
for(n=0; n<Ml; n++){ 

unsigned char* x = malloc(16); 
for(m=0; m<16; m++) x[m] = a[n*16+m]; 

c [n] = *(double complex*)x; 

} 

double complex *cc = calloc(M2, sizeof(double complex)) 
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for(n=0; n<M2; n++){ 

unsigned char* x = malloc(16); 

for(in = 0; in<16; m + + ) x [m] = aa[n*16 + m]; 

cc[n] = *(double complex*)x; 

} 

double complex *uu = calloc(NN, sizeof(double complex)) 
double complex *vv = calloc(NN, sizeof(double complex)) 
double complex *r = calloc(NN, sizeof(double complex)); 
for(m=0; m<Ml; m++){ 

double norm = fnorm(NN, pu) ; 

for(n=0; n<NN; n++) uu[n] += c[m]*pu[n]/norm; 

cprod(N, pu, u, r); 

for(n = 0; n<NN; n + + ) pu[n] = r[n] ; 

} 

for(m=0; m<M2; m++){ 

double norm = fnorm(NN, pv) ; 

for(n = 0; n<NN; n + + ) vv [n] += cc[m]*pv[n]/norm; 

cprod(N, pv, V, r); 

for(n = 0; n<NN; n + + ) pv[n] = r[n] ; 

} 

double complex *b = calloc(NN, sizeof(double complex)); 
for(n=0; n<NN; n++){ 

unsigned char* x = malloc(16); 
for(m=0; m<16; m++) x[m] = B[n*16+m]; 
b [n] = *(double complex*)x; 

} 

cprod(N, uu, b, pu); 
cprod(N, pu, vv, r); 
for(n=0; n<NN; n++){ 

double X = significand(creal(r[n])); 
if (x > 0) X = X - 1; else x = x + 2; 
unsigned fr = (unsigned) (x*4294967295) ; 
double y = significand(cimag(r[n])); 
if(y > 0) y = y - 1; else y = y + 2; 
unsigned fi = (unsigned)(y*4294967295); 
unsigned char* rr = (unsigned char*)&fr; 
unsigned char* ri = (unsigned char*)&fi; 
for(m=0; m<4; m++) S[n*8+m] = rr[m]; 

for(m = 0; m<4; m + + ) S [n*8+m+4] = ri [m] ; 

printf ( "7«u\t7«u\n" , f r , fi); 

} 

free(u); free(v); free(c); free(cc); 
free(pu); free(pv); free(uu); free(vv); 
free(b); free(r); 

} 

void public_key(unsigned Ml, unsigned M2, unsigned N, 
unsigned char *a, unsigned char *aa, 
unsigned char *U, unsigned char *V, 
unsigned char *A){ 
unsigned m, n, i, j, NN = N*N; 

double complex *u = calloc(NN, sizeof(double complex)); 
double complex *v = calloc(NN, sizeof(double complex)); 
double complex *pu = malloc(NN*sizeof(double complex)); 
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double complex *pv = malloc(NN*sizeof(double complex)); 
if (U[0] >= V[0]){ 

for(n=0; n<NN; n++){ 

unsigned char* x = malloc(16); 
for(m=0; m<16; m++) x[m] = U[n*16+m]; 
pu[n] = u[n] = *(double complex*)x; 

} 

for(n=0; n<NN; n++){ 

unsigned char* x = malloc(16); 
for(m=0; m<16; m++) x[m] = V[n*16+m]; 
pv[n] = V[n] = *(double complex*)x; 

} 

} 

else{ 

for(n=0; n<NN; n++){ 

unsigned char* x = malloc(16); 
for(m=0; m<16; m++) x[m] = V[n*16+m]; 
pu[n] = u[n] = *(double complex*)x; 

} 

for(n=0; n<NN; n++){ 

unsigned char* x = malloc(16); 
for(m=0; m<16; m++) x[m] = U[n*16+m]; 
pv[n] = V[n] = *(double complex*)x; 

} 

} 

double complex *c = calloc(Ml, sizeof(double complex)); 
for(n=0; n<Ml; n++){ 

unsigned char* x = malloc(16); 
for(m=0; m<16; m++) x[m] = a[n*16+m]; 

c[n] = *(double complex*)x; 

} 

double complex *cc = calloc(M2, sizeof(double complex)) 
for(n=0; n<M2; n++){ 

unsigned char* x = malloc(16); 
for(m=0; m<16; m++) x[m] = aa[n*16+m]; 

cc[n] = *(double complex*)x; 

} 

double complex *uu = calloc(NN, sizeof(double complex)) 
double complex *vv = calloc(NN, sizeof(double complex)) 
double complex *r = calloc(NN, sizeof(double complex)); 
for(m=0; m<Ml; m++){ 

double norm = fnorm(NN, pu); 

for(n=0; n<NN; n++) uu[n] += c[m]*pu[n]/norm; 

cprod(N, pu, u, r); 

for(n = 0; n<NN; n + + ) pu [n] = r[n] ; 

} 

for(m=0; m<M2; m++){ 

double norm = fnorm(NN, pv); 

for(n = 0; n<NN; n + + ) vv[n] += cc [m]*pv[n]/norm; 

cprod(N, pv, V, r); 

for(n = 0; n<NN; n + + ) pv [n] = r[n] ; 

} 

cprod(N, uu, vv, r); 



for(n=0; n<NN; n++){ 

double complex cr = r [n] ; 

unsigned char* rr = (unsigned char*)&cr; 

for(m=0; m<16; m++) A[n*16+m] = rr[m]; 

pr int f ( " 7o+. 161f\t 7o+. 16 If \n " , creal(cr), cimag(cr)); 

} 

free(u); free(v); 
free(c); free(cc); 
free(pu); free(pv); 
free(uu); free(vv); free(r); 

} 

void matrix_public_key(unsigned N, unsigned char *U){ 
unsigned n, m; 
for(n=0; n<N/16; n++){ 

double complex c = drand()*2 - 1 + l*(drand()*2 - 1); 
unsigned char* array = (unsigned char*)&c; 
for(m=0; m<16; m++) U[n*16+m] = array[m]; 
pr int f ( " 7«+. 161f\t 7«+. 16 If \n " , creal(c), cimag(c)); 

} 

} 

void private_key(unsigned N, unsigned char *a){ 
unsigned m, n; 
for(n=0; n<N/16; n++){ 

double complex c = drand()*2 - 1 + l*(drand()*2 - 1); 
unsigned char* array = (unsigned char*)&c; 
for(m=0; m<16; m++) a[n*16+m] = array[m]; 

pr int f ( " 7«+. 161f\t 7«+. 16 If \n " , creal(c), cimag(c)); 

} 

} 

void sprint(unsigned N, unsigned char *x){ 
unsigned n; 
for(n=0; n<N; n++){ 
if(n 7. 32 == 0) 

printf ( "\n7«02x" , x [n] ) ; 
else 

pr intf ( " 7«02x " , x [n] ) ; 

} 

printf("\n"); 

} 

unsigned char check_key(unsigned N, unsigned char *x, unsigned char *y){ 
unsigned n, c = 0; 
for(n=0; n<N; n++){ 
if(x [n] ! = y[n] ) { 

c = 1; 
break; 

} 

> 

return c; 


9 



mainCint argc , char *argv[]){ 

// Dimension of matrices, NxN 
unsigned N = 4, NN = N*N; 

printf("Dimension of matrices\n"); 

printf ("-\n") ; 

printf ("N = 7.u\tNN = 7.u\n", N, NN); 

// Length of the secret key 
unsigned K = NN*8; 

printf("\nLength of the secret key\n"); 

printf ("-\n"); 

printf ("K = 7ou\n", K) ; 

// Initialize the random number generator 
srand ( (uns igned char) time(O) + getpidO); 

// Alice: private key vectors 

unsigned Ml = NN + rand()7.NN, M2 = NN + rand()7.NN; 
unsigned char *a = malloc(Ml*16); 
unsigned char *aa = malloc(M2 * 16) ; 
printf("\nAlice: private key a\n"); 

printf ("-\n") ; 

printf ("\nMl = 7«d complex numbers\n", Ml); 
private_key(Ml*16, a); 

pr intf ("\nHex representation of a, 7«u bytes", Ml*16); 
sprint(Ml * 16 , a); 

printf("\nAlice: private key aa\n"); 


printf ("-\n") ; 

printf ("\nM2 = 7«d complex numbers\n", M2); 


private_key(M2*16, aa); 

pr intf ("\nHex representation of aa , 7«u bytes", M2*16) 
sprint(M2*16 , aa); 

// Bob: private key vectors 


unsigned J1 = NN + rand()7.NN, J2 = NN + rand()7.NN; 
unsigned char *b = malloc(J1* 16) ; 
unsigned char *bb = malloc(J2* 16) ; 
printf("\nBob: private key b\n"); 

printf ("-\n") ; 

printf ("\nJl = 7«d complex numbers\n", Jl); 


private_key(J1* 16 , b); 

pr intf ("\nHex representation of b, 7«u bytes", Jl*16); 
spr int ( Jl * 16 , b) ; 

printf("\nBob: private key bb\n"); 


printf ("-\n") ; 

printf ("\nJ2 = 7«d complex numbers\n", J2); 


private_key(J2*16 , bb); 

pr intf ("\nHex representation of bb , 7«u bytes", J2*16) 
sprint(J2*16 , bb); 
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// Alice: matrix public key 

unsigned char *U = malloc(NN* 16) ; 

printf("\nAlice: matrix public key U\n"); 

printf ("-\n"); 

printf ("\nN*N = 7,d complex numbers\n", NN); 
matrix_public_key(NN*16, U); 

pr intf ("\nHex representation of U, 7oU bytes", NN*16); 
sprint(NN*16, U); 

// Bob: matrix public key 

unsigned char *V = malloc(NN* 16) ; 

printf("\nBob: matrix public key V\n"); 

printf ("-\n"); 

pr intf ("\nN*N = 7od complex numbers\n", NN); 
matrix_public_key(NN*16, V); 

pr intf ("\nHex representation of V, 7«u bytes", NN*16); 
sprint(NN*16 , V) ; 

// Alice: public key 

unsigned char *A = malloc(NN* 16) ; 

printf("\nAlice: public key A\n"); 

printf ("-\n") ; 

pr intf ("\nN*N = 7«d complex numbers\n", NN); 
public_key(Ml, M2, N, a, aa, U, V, A); 

pr intf ("\nHex representation of A, 7«u bytes", NN*16); 
sprint(NN*16 , A) ; 

// Bob: public key 

unsigned char *B = malloc(NN* 16) ; 
printf("\nBob: public key B\n"); 

printf ("-\n") ; 

pr intf ("\nN*N = 7od complex numbers\n", NN); 
public_key(J1, J2, N, b, bb, U, V, B); 

pr intf ("\nHex representation of B, 7«u bytes", NN*16); 
sprint(NN*16 , B) ; 


// Alice: secret key 
unsigned char *Sa = malloc(K); 
printf("\nAlice: secret key Sa\n"); 
printf ("-\n") ; 

pr intf ("\n2 *N*N = 7od unsigned integers on 4 bytes\n", 2*NN) 
secret_key(Ml, M2, N, K, a, aa, U, V, B, Sa); 


pr intf ("\nHex representation of Sa , 7«u bytes", K) ; 
sprint(K , Sa) ; 

// Bob: secret key 
unsigned char *Sb = malloc(K); 
printf("\nBob: secret key Sb\n"); 
printf ("-\n") ; 

pr intf ("\n2 *N*N = 7od unsigned integers on 4 bytes\n", 2*NN) 
secret_key(J1, J2, N, K, b, bb, U, V, A, Sb); 
pr intf ("\nHex representation of Sb , 7«u bytes", K) ; 
sprint(K , Sb) ; 
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// Check the secret keys 

printf("\nCheck the secret keys\n"); 

printf ("-\n") ; 

if(check_key(K, Sa, Sb) == 0){ 
printf("Result; Sa = Sb\n"); 

} 

else{ 

printf("Result: Sa != Sb\n"); 

} 

return 0; 


Appendix 2 

For illustration purposes, the results obtained for one instance run are given below. 
Size of matrices 


N = 4 NN = 16 

Length of the secret key 


K = 128 

Alice; private key a 


Ml = 24 complex numbers 
+0.8689850983830614 
-0.1840404501840445 
-0.5624424497754633 
+0.2681173127751544 
+0.6868307886138416 
-0.5409692556398616 
+0.0340424836674749 
+0.4317577723105932 
+0.0672014622485229 
+0.9595978984392544 
+0.7195946921855143 
+0.5863311031238485 
+0.7634357203607813 
+0.2165438627036125 
+0.2439741527589443 
+0.5183932032859373 
-0.4087175890232525 
-0.5239555517482943 
-0.0022945854777570 
+0.4798449924309158 
-0.4060460210598060 
+0.6693677004712948 
-0.1540452316059945 
+0.3646444266501574 


-0.8548606934416402 
-0.9778283575338482 
-0.3926427336457170 
-0.2169790138199983 
-0.9872768015372509 
-0.8628458723853144 
-0.3506396088873717 
+0.7887667341912554 
+0.2782204334161911 
-0.0538938046045783 
+0.8729643758373127 
-0.7404569811244905 
+0.6367243534294711 
+0.5814044504531037 
+0.7957204292763893 
-0.9873300964591410 
+0.9764747939005469 
-0.5312443242637066 
-0.0519642428515297 
+0.1833511453978887 
+0.8571984084524591 
+0.4657471765812942 
-0.9389158129164289 
+0.5890737755058091 
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Hex representation of a, 384 bytes 

ea8748d6b9ceeb3f9df71ed0045bebbf483e5731a38ec7bf5eb216b25e4aefbf 
Ibbef54e87ffelbf204e01fd0e21d9bf90157784d528dl3f903a23e4f7c5cbbf 
526fdf8f84fae53f9576d684c597efbf5faba3cl9e4felbfcl456cf26e9cebbf 
c0860d34046eal3f9e9f361del70d6bf3c8af759ebaldb3ffac5f0bb933de93f 
00cd99721d34bl3fc07ea6135dcedl3ff26be3a606b5ee3fa0a47700f897abbf 
aafba972eb06e73ffe6a99fc52efeb3fd06cll7239c3e23f54c9c0d6d2ble7bf 
3a8771bf106ee83f7c7b51c00b60e43f68133b94b5b7cb3fe43c8e81dd9ae23f 
f0a095878b3acf3f9c0990b08a76e93ff29dd257ad96e03fcfc054493598efbf 
dee889dl6d28dabf627a2511483fef3f4221ea6e3ec4e0bf60bedcl8f4ffeObf 
0008b73218cc62bf60aa37a70e9baabf3cf868c5c7b5de3f688395e20c78c73f 
dcdb4773a8fcd9bfea914f5b2b6eeb3fe2bbd0cf756be53f64bee73ecdcedd3f 
741eed0fclb7c3bfcf8bc52c990beebfe463c8935556d73f7a51173fbld9e23f 

Alice; private key aa 


M2 = 22 complex numbers 

-0.5118704412976753 -0.0500573637169195 
+0.7604183635929800 -0.4931928538884157 
-0.3375783216472882 +0.1598131473912503 
+0.1632397906851764 +0.3421602938986694 
-0.3620384387433089 +0.3855473826440021 
-0.8872378358483427 +0.5676916094065314 
+0.7025298861832101 +0.2206711130955801 
-0.6247065410609323 +0.3230649032942832 
+0.6917058497340922 +0.0231011495289499 
-0.0152488020307338 +0.8275970245954227 
+ 0.3305744647758120 -0.78 97 5 200 94890274 
-0.2482330123079071 +0.7894395461742834 
-0.9737185128877407 -0.7489919758449641 
+0.4474754300516699 +0.6970578097454476 
+0.7274738114424169 +0.0977605819469343 
-0.3061669770440641 -0.7828203718630750 
-0.2298109276249133 -0.0338097119771229 
+0.7811974836736504 +0.8376630850857776 
+0.6003193600063923 +0.4021131457350087 
-0.6334603792724227 +0.6055862039104507 
+0.6752938616270538 +0.7418148947846244 
-0.6433983716975780 -0.4681352436902144 

Hex representation of aa, 352 bytes 

0c34a51e3e61e0bf80212d68leala9bf261c5de45855e83ff29584c27890dfbf 
3213d41ae29ad5bf6814c2d8cl74c43f3009339d0ae5c43fcc35124af4e5d53f 
c8089345a32bd7bfb8ebe0edceacd83f23c2179a4064ecbf48abl398872ae23f 
dcd2b3f41f7be63f786ef576f33ecc3f0c876e9298fde3bf2c99886al8add43f 
04ea614e7422e63f404ae8e6d3a7973f80913390c33a8fbf9aef5cclac7bea3f 
50b5c6cc2128d53f68949601a645e9bfl034d36el9c6cfbfa89alfb91643e93f 
6d9a0bbab328efbfb00e2705bef7e7bf247875fc6fa3dc3f18e6082e4c4ee63f 
ceee9a287747e73f80cd4c66d606b93f7a3e61603d98d3bf436cf94edd0ce9bf 
5cca34c9716acdbf503575ae814falbffal083dd91ffe83ff05070d022ceea3f 
66424cf2d035e33f8c4e8ec638bcd93f902d89b34e45e4bf96889651f660e33f 
7el65cdf019ce53ffc231997f2bce73f38b4972eb896e4bf18497086edf5ddbf 
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Bob: private key b 


J1 = 31 complex numbers 
-0.8131572898736764 
+0.1298227154186107 
-0.7446356426328338 
-0.1602102571093867 
-0.0321974180566288 
+0.6083909207849130 
+0.1000333272410052 
-0.4108327516259943 
-0.4353843212160295 
-0.9065787271243225 
+0.6038810599955362 
-0.7069626958783282 
-0.9291368614924062 
+0.1682809573220776 
-0.3759606473197272 
+0.4097225812747536 
+0.0968815767756941 
+0.7941227780332327 
-0.9919780830175842 
-0.5993550677175068 
+0.9217383104992827 
+0.6513040121506304 
+0.5596626098265201 
-0.4966598369750350 
+0.3960260341933464 
-0.0567990502670347 
+0.8749646723409310 
+0.1605486771157856 
-0.6534623397058436 
-0.1318483317593858 
-0.4699044210083154 


-0.9495491783647080 
+0.7302350146477263 
+0.9468203098361780 
+0.0153968761011387 
-0.1509698059198559 
-0.6980298260647722 
-0.8515098187419327 
+0.7698810134887732 
+0.8627621712051379 
+0.9290429557182074 
-0.3380758062243528 
-0.9937368694489492 
-0.1109532512441092 
+0.9704573627930386 
+0.8195097335932378 
+0.6411222197176920 
+0.9892165586232851 
+0.6687766533299828 
+0.5390884771002475 
+0.8416746390522734 
+0.5618991062378083 
-0.8492378708487287 
+0.3583335523844433 
+0.4885499630912138 
-0.0562765497713249 
-0.2094754124106146 
-0.4118473497446520 
+0.2522820859568238 
-0.2528823245864887 
-0.8276937454948939 
-0.7268149655182780 


Hex representation of b, 496 bytes 
a75dd06f6205eabf97a460f5b462eebfc81b80de079ec03f 
286891200ed4e7bf4aeb3dlb5a4cee3f6c60620bc581c4bf 
b04123132e7ca0bfd8fb8d85fa52c3bf7e8b7e3cf077e33f 
609b02bdc89bb93f0d06f784913febbf5eff1674154bdabf 
acl8343356dddbbf26aabl69bf9beb3fdbf10764b102edbf 
3a596b5ffe52e33ffe5ad3b408a3d5bfe042493b709fe6bf 
61c6333a7dbbedbf30737aa96e67bcbf58731efc3a8ac53f 
d88e343fbd0fd8bf102cll7a6c39ea3f487fcOOfe538da3f 
2026d6233bcdb83f8e5bfe7ba9a7ef3f342al52c7469e93f 
f6181dd248beefbfe678597a3640el3fdc3bdladea2de3bf 
20126257el7eed3f302d9ed513fbel3ff018fe827bd7e43f 
16f7bf8fcle8el3f743e23daefeed63f8cd8425746c9dfbf 
68664e947d58d93f805f99a947d0acbff079aedec314adbf 
4a829be9b5ffeb3f543e86fcb45bdabf306ed0eadb8cc43f 
238946da29e9e4bf9e777558392fd0bffc2778f867e0c0bf 
787flefee912debfbl9e64751142e7bf 


4acb49d2155ee73f 
000785ba65888f3f 
a09552a54256e6bf 
16dcd781dda2e83f 
ae71ec4ab8baed3f 
e69a6343blccefbf 
e6786b99fc0def3f 
2ea5cdbel284e43f 
ec2dcc4b9e66e53f 
8e4413a7ffeeea3f 
203d3ae6f42cebbf 
d0147cl06744df3f 
6852cfIel7d0cabf 
2c4723c36325d03f 
5d230098777ceabf 
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Bob: private key bb 


J2 = 17 complex numbers 
+0.2395015183561122 
-0.0465769789069043 
+0.7537951322336816 
-0.1563436965923156 
-0.0856898852384244 
+0.5493668891763306 
-0.2105619723294515 
+0.2936094795551767 
-0.4165447796757208 
-0.5891911056850502 
+0.4081439302681709 
-0.6185220295108234 
-0.6507796759984934 
-0.3987045968811376 
+0.3471997888579212 
+0.4039709032807668 
-0.9501367908238625 


+0.3971994682370852 
-0.3663425200453365 
+0.8256763657073602 
-0.1521531633898803 
+0.9043694082667142 
+0.0864919672072155 
+0.2586894465691962 
-0.1546671663785663 
-0.6470306531837624 
+0.1686385513920701 
+0.4335022113121134 
+0.7933638301919215 
-0.7439947848693523 
-0.4433641387621932 
-0.6838942312678632 
-0.2897329644110010 
+0.1739202672557096 


Hex representation of bb, 272 bytes 


b846575afca7ce3fa0478451b76bd93f 
2c7claf8161fe83f4a6279d7f06bea3f 
4870b2b6c5efb5bf6444001d98f0ec3f 
14d777d8blf3cabf4071352e5e8ed03f 
36b6816fabaSdabf82eldda079b4e4bf 
400324b8071fda3fc011150f80bedb3f 
86142ae62fd3e4bfcfllad26cecee7bf 
04aa94768538d63f39fea62776e2e5bf 
84136a458567eebff8c7fcfl0443c63f 


30a34dl2f0d8a7bf08a7aee52772d7bf 
f45de6fbll03c4bf98095f3ecl79c3bf 
2690d0de6994el3fd08d856a5624b63f 
4c0blf6a7fcad23f0c33ae3a22ccc3bf 
255d404ea7dae2bfd87089b3f295c53f 
565813b6eecae3bf7eldl08b3c63e93f 
6ca317496084d9bf18fb0cfbl360dcbf 
201888c6a8dad93ff6f9ad21fc8ad2bf 


Alice: matrix public key U 


N*N 

+0.3368737329532423 
+0.2391556223941222 
-0.8836268003784012 
-0.4046797065655433 
+0.2883828453162156 
+0.8417251050926819 
-0.6391474565256796 
-0.2071290143032435 
+0.7152330705414824 
+0.9153039464129042 
+0.1244033889339511 
-0.0354984985378332 
+0.5417379064036945 
-0.1986953363917494 
-0.3848289466173010 
-0.3802923558073766 


+0.5243043640910983 
-0.3747238699602812 
-0.0031821031027317 
-0.0835783218659490 
-0.6358192650617027 
-0.4126508745199150 
+0.2763919025872350 
+0.2331221310471723 
-0.1286121438584548 
-0.1252376190599107 
-0.4249777151124192 
-0.5027293312737142 
-0.0836606229409183 
-0.1115330181161580 
-0.3708409810614971 
+0.5728729606643514 


16 complex numbers 
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Hex representation of U, 256 bytes 
68997ad8568fd53fa67bldf219c7e03f002c6bc4a69cce3f 
bfdl2fb6ab46ecbf008a47985all6abfda7243b645e6d9bf 
a0135752dd74d23fa66ab3a4al58e4bf240cd37c69efea3f 
162ee35de573e4bff8729fa967b0dl3f4cl93elb3483cabf 
58c9df7630e3e63f0cc9dedb5c76c0bff8c777802b4aed3f 
909efb86e6d8bf3fOOf2fbbad532dbbfbOdf8df4db2ca2bf 
6e3ce0bbea55el3ff0857e57c86ab5bfd0cbd549d96ec9bf 
c6a0119709ald8bf6e3e6bcfdbbbd7bf7622c7bfb556d8bf 


42a3a0d379fbd7bf 
b89dll8f6365b5bf 
3aa5ae36df68dabf 
a824692cf2d6cd3f 
30ad0b4bc907c0bf 
Ocf091d25bl6e0bf 
c0aa3b896d8dbcbf 
5clddaacf954e23f 


Bob: matrix public key V 


N*N 

+0.9409248020558381 
-0.7158705844244213 
+0.3557956546476002 
+0.2466279650314691 
+0.5482843655448633 
+0.6265073683117459 
-0.3367309500069268 
+0.5288800137059331 
+0.6000563534476080 
-0.2284014936471640 
+0.5991160201006454 
-0.5344290793776173 
-0.5174786942983332 
-0.3930733586311385 
-0.0401902524752414 
+0.1879347693971776 


+0.5377867567981318 
+0.2893294067383878 
+0.8879058945032454 
+0.1282564850152923 
-0.4268184961805279 
-0.8281711179997823 
+0.5117462362356848 
+0.6497204598425779 
+0.6179174079884797 
+0.3260335651618362 
+0.2927185726013939 
-0.5523086828897967 
+0.6341098574986008 
-0.6365556865422852 
-0.6492929052234647 
-0.4468276153836913 


16 complex numbers 


Hex representation of V, 256 bytes 
34679a540elcee3fac6f95928c35el3f76ae886d69e8e6bf 
784d31235bc5d63f66255a9fb969ec3fa06a0b528191cf3f 
a68b5da78b8bel3fd4169b86fe50dbbf6243332e590ce43f 
322b75f8ff8cd5bfa8778fa43960e03fIe004cc795ece03f 
5410ba61a933e33f48115ebafac5e33f443bc998423ccdbf 
42f1la5cf52be33f744el0aee6bbd23f05aa3e030blaelbf 
806b8c7a2f8feObf02e380cla04ae43fc2f90f291d28d9bf 
f08fff17dl93a4bfIdb52eea01c7e4bf88cc2blc3f0ec83f 


e898ed7c5f84d23f 
OOff5160b56ac03f 
55a069b76080eabf 
4alcd28f82cae43f 
002dfIe2bbddd43f 
54df494283acelbf 
ac05f907aa5ee4bf 
9e72cldad298dcbf 


Alice: public key A 


N*N 

-0.2856153551337328 
+0.0053028357754545 
-0.6848378625559008 
-0.1177416781224439 
-0.1356525197048316 
-0.0119241534869403 
-0.1496198348578027 
+0.1844067261375284 


-0.0836362385140290 
-0.6553666223075714 
+1.0033521606714080 
+0.1881147013745262 
-0.4999190536353746 
+0.1657257251876224 
-0.2369035535027847 
+0.0291454067840425 


16 complex numbers 
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+0.1145145558807253 
-0.1005465316239941 
-0.4039634042904244 
-0.2481990593712110 
+0.1218629445687005 
+0.5246569994564118 
-0.5467684445921881 
-0.2094934354311639 


+0.0391101989223119 
-0.6368000194347039 
+1.0076922004908939 
+0.2515002673868378 
+0.0144359211753338 
+0.0595147923884747 
-0.4633383063630160 
-0.2169002141819963 


Hex representation of A, 256 bytes 
4c3a62a08547d2bf309f2d3d2f69b5bfa0c7b4236db8753f 
34afd71731eae5bf2f46c7feba0df03f747ce9905124bebf 
a249e0cfOf5dclbf90b7807cacfedfbf803612cab06b88bf 
Ic0ac624be26c3bffeel4c0bdb52cebf9128f6bca39ac73f 
42766 c7 0d3 5 0bd3fc26d9f8 2400 6a43ff8200del6abdb9bf 
22538d5289dad9bfclae58db8Ilff03f2adb729dfcc4cfbf 
f34e84f16832bf3f857faf8a94908d3f850cc979fdc9e03f 
78438089207felbf2e0b34b655a7ddbfl5a7f34eaed0cabf 


9e6e366cc3f8e4bf 
6273267d2414c83f 
4cb3e4248036c53f 
c843dd234bd89d3f 
120e326faa60e4bf 
7b26d7939418d03f 
e09adf0db978ae3f 
dbla34df62c3cbbf 


Bob: public key B 


N*N 

-2.4633752471535728 
-1.7414343240690795 
-0.7089011929436275 
+0.2888098643832883 
-0.4447799179419522 
-0.1946598847601872 
-0.0495878578361504 
+0.0332765292089976 
-2.0814149405134770 
-1.6706358918192217 
+0.2365166391262642 
+0.6503887269428416 
+0.6959317836726924 
+0.7908533006696801 
-0.6917520195247205 
-0.5250685199630863 


-0.1602678155173908 
-1.0301803412988897 
+2.6033653995146535 
+1.2726176834487517 
+0.2067805379384524 
-0.0414232823029533 
+0.4089135120393886 
+0.0188205856998400 
+0.5097364399723628 
-0.2970409670111223 
+2.1543529375391444 
+0.8676760384430399 
-0.6902765198764746 
-0.2909750558709320 
-0.9240611824254524 
-0.1903208570713778 


16 complex numbers 


Hex representation of B, 256 bytes 
7367e214feb403c0a6d01fela783c4bfl223e03ceadcfbbf 
7a39f98d51afe6bf96a0133dbld304407275925edc7bd23f 
9e9b5e304677dcbf678df2dfc877ca3f57ae71779deac8bf 
80b35e679463a9bfde814694a32bda3f20a45a099f09al3f 
134a57e0bca600c0315868cbc24fe03f97326eb3ecbafabf 
ede9005f2d46ce3f602d96641d3c0140bl02fd04fccfe43f 
8ce063bbl245e63fecldc2c8bel6e6bff4eac994ab4ee93f 
bf9e9921d522e6bf2fa7c0cle891edbf8bdl2c7f5ccde0bf 


0ac6ad619e7bf0bf 
26942b5ca45cf43f 
707d91b56e35a5bf 
821d4a20b445933f 
99a3b81db802d3bf 
9d57148a00c4eb3f 
Ibb73ad7559fd2bf 
13236f106f5cc8bf 


Alice: secret key Sa 


2*N*N = 32 unsigned integers on 4 bytes 
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3532271144 

872140306 

3218067163 

1175145558 

2871309341 

380671065 

1259174687 

3387166373 

1308192945 

212774227 

4247529633 

4198654425 

2079979731 

3773268400 

1322895297 

4146361598 


1241794686 

1131329665 

1418852160 

3891638795 

18700978 

746843244 

3031792119 

3348009390 

820672386 

4220132680 

3762180882 

3427153554 

1375907524 

1991239246 

3617498003 

273775259 


Hex representation of Sa, 128 bytes 

282e8ad27e48044al2cefb3381b86e43dbcecfbf40f79154564c0b460bb2f5e7 
Idb424abb25ald015994b0166cec832clf7b0d4bf779b5b4a50ee4c9ae918ec7 
bl70f94d8277ea3053adae0c481d8afbal282cfdl2533ee0d96142fa923646cc 
d3f8f97bc4ae0252b081e7e04ee6af76clc7d94e93a39ed7fe7424f79b7a5110 

Bob: secret key Sb 


unsigned integers 


2*N*N = 32 

3532271144 

872140306 

3218067163 

1175145558 

2871309341 

380671065 

1259174687 

3387166373 

1308192945 

212774227 

4247529633 

4198654425 

2079979731 

3773268400 

1322895297 

4146361598 


1241794686 

1131329665 

1418852160 

3891638795 

18700978 

746843244 

3031792119 

3348009390 

820672386 

4220132680 

3762180882 

3427153554 

1375907524 

1991239246 

3617498003 

273775259 


on 4 bytes 


Hex representation of Sb, 128 bytes 

282e8ad27e48044al2cefb3381b86e43dbcecfbf40f79154564c0b460bb2f5e7 
Idb424abb25ald015994b0166cec832clf7b0d4bf779b5b4a50ee4c9ae918ec7 
bl70f94d8277ea3053adae0c481d8afbal282cfdl2533ee0d96142fa923646cc 
d3f8f97bc4ae0252b081e7e04ee6af76clc7d94e93a39ed7fe7424f79b7a5110 


Check the secret keys 


Result: Sa = Sb 
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